BOUNDARY CHARACTERISTIC POINT REGULARITY FOR 

NAVIER-STOKES EQUATIONS: BLOW-UP SCALING AND 

PETROVSKII-TYPE CRITERION (A FORMAL APPROACH) 

V.A. GALAKTIONOV AND V. MAZ'YA 

^— ^ ' Abstract. The three-dimensional (3D) Navier-Stokes equations 

(0.1) Ut + (u- V)u = -Vp + Au, divu = in Qq, 

where u — [u, v, w]'^ is the vector field and p is the pressure, are considered. Here, 
Qo CM.^ X [—1,0) is a smooth domain of a typical backward paraboloid shape, with the 
vertex (0, 0) being its only characteristic point: the plane {t = 0} is tangent to dQo at 
the origin, and other characteristics for t e [0,-1) intersect dQo transversely. Dirichlet 
Ph ' boundary conditions on the lateral boundary dQo and smooth initial data are prescribed: 

(0.2) u = on dQo, and u(a;, -1) = Uo(a;) in Qo n {t = -1} (divuo = 0). 

Existence, uniqueness, and regularity studies of (jO.ip in non- cylindrical domains were 
initiated in the 1960s in pioneering works by J.L. Lions, Sather, Ladyzhenskaya, and 
Fujita-Sauer. However, the problem of a characteristic vertex regularity remained open. 
In this paper, the classic problem of regularity (in Wiener's sense) of the vertex (0, 0) 
for (|0.1|) . (|0.2p is considered. Petrovskii's famous "2v^IogTog-criterion" of boundary 
A^ • regularity for the heat equation (1934) is shown to apply. Namely, after a blow-up scaling 

^\ , and a special matching with a boundary layer near dQo , the regularity problem reduces 

r^ ' to a 3D perturbed nonlinear dynamical system for the first Fourier-type coefficients of the 

solutions expanded using solenoidal Hermite polynomials. Finally, this confirms that the 
nonlinear convection term gets an exponentially decaying factor and is then negligible. 
Therefore, the regularity of the vertex is entirely dependent on the linear terms and 
hence remains the same for Stokes' and purely parabolic problems. 

Well-posed Burnett equations with the minus bi-Laplacian in (jO.ip are also discussed. 
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1. Introduction: vertex regularity for the Navier-Stokes equations 

1.1. Navier— Stokes equations inside a non-cylindrical backward paraboloid: 
first history since 1960s. We consider 3D Navier-Stokes equations (the 3D NSEs) 

(1.1) ut -h (u • V)u = -Vp + Au, divu = in Qo, 
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where u = [u, v, w]'^{x, t) is the vector field and p = p{x, t) is the corresponding pressure. 
The NSEs f ll.ip are posed in a smooth non- cylindrical domain 

QoCR-^x [-1,0) 

of a typical backward paraboloid shape, with the vertex (0, 0) being its only characteristic 
point the plane {t = 0} is tangent to OQq at the origin. No characteristic points of 
dQo are assumed to exist for t e [—1,0), i.e., other characteristics {t = r}, for any 
r G [—1,0), intersect dQo transversely, in a natural sense. Next, the zero Dirichlet 
boundary conditions on the lateral boundary dQo and smooth initial data at t = — 1 are 
prescribed: 

u = on dQo, and 
(1.2) 

u(x, — 1) = uo(a;) in Qo H {t = — 1}, where divuo = 0. 

The questions of solvability, uniqueness, and regularity for the Navier-Stokes equations 
in non- cylindrical (and non-characteristic) domains, i.e., in our case, up to the vertex, 
for t < —6o < 0, were actively studied since the 1960s. J.L. Lions began this study in 
1963; see references in his classic monograph jHJ Ch. 3] concerning elliptic regularization- 
penalization methods; as well as Fujita-Sauer [151 [16] (1969) as one of the first such study 
of weak solutions via a penalization. Another alternative, as was pointed out in jUJ Ch. 3, 
§ 8.1], is a "rather careful using" Galerkin methods with time dependent basis functions; 
see Sather [64j (1963). In 1968, Ladyzhenskaya [39] proved local existence (global for 
N = 2 and for small initial data if A^ = 3) and uniqueness of strong solutions for time- 
dependent domains using a different method. See Neustupa [SH] for more recent results, 
references, and other related problems. 

However, the problem of regularity of a characteristic boundary point for the NSEs in 
any dimension N > 2 was not addressed elsewhere and remained open. Naturally, in 
order to proceed with regularity issues concerning the paraboloid vertex (0,0), we have 
to assume that a unique smooth bounded solution of (II. ip . (II. 2p exists in Qq, i.e., with 
no L°°-blow-up for t < Qj. In particular, as is well-known (see [2S1 SSI ES] ) , global smooth 
solutions always exist for sufficiently small initial data, so we can directly proceed with 
the vertex regularity, at least, for this class of solutions. 

1.2. Regularity of the characteristic paraboloid vertex. Thus, the classic problem 
of regularity (in Wiener's sense, see [47j) of a boundary characteristic point for the NSEs 
problem (II. ip . (II. 2p is under consideration. 

Definition (vertex regularity/irregularity). According to Wiener [70], the vertex 
(0,0) of the given backward paraboloid Qo for the NSEs problem (II. ip . (I1.2p is regular 
if, for any bounded data Uo(a;), 

(1.3) u(0,0") = 0, 

and irregular otherwise, i.e., at least for some initial data, (II. 3p fails. 



But the solution is formally allowed to blow-up at the vertex (0, ). 
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The boundary and other regularity issues for the Navier-Stokes equations in M? and M^ 
have been and remain key and very popular in modern mathematical literature, since 
J. Leray's seminal papers in 1933-34 ^21 HS]. Among various regularity and partial 
regularity results for the NSEs, the boundary regularity properties in piecewise smooth 
or Lipschitz domains and those with thin channels, or other non-regular domains (as we 
will show, such settings are key for our study) always played a special role. Mentioning 
Kondratiev's first study of 1967 [33], we refer to advanced results, further references, and 
reviews in recent papers [TJ |3T1 |38l SU EOj |52l [56] and [71] . See also [STJ [57] for the related 
linear Stokes problem 

(1.4) ut = — Vjo + Au, divu = in Qo, u(0,a:) = uo(x) (divuo = 0). 

Concerning compressible flows and other related problems, see a good survey in |38j . 
where 2D NSEs in a polygon domain with a convex vertex were studied. 

Note that, and this is key for us in what follows, J. Leray in [121 113] actually posed 
a deep problem on both backward and forward continuation phenomena, which sound 
modern and advanced nowadays for general nonlinear PDE theory: 

Leray's blow-up scenario: self-similar blow-up as t — ?■ T^ [t < T) 
and similarity collapse of this singularity as t — >■ T+ (t > T) ; 

see his precise statements and a discussion on these principal issues in [18i § 2.2]. In this 
connection, such "backward blow-up scaling approaches" will be key later on. 

According to our approach, we deal with a typical asymptotic problem of clarifying a 
generic behaviour of solutions near a "blow-up singularity" (0,0) (a "micro-scale struc- 
ture" of nonlinear PDEs involved). Of course, the vertex regularity problem setting 
essentially and crucially depends on the a priori given shape of the prescribed backward 
paraboloids, which affects our methods of matched blow-up expansions. Anyway, we hope 
that our blow-up analysis on shrinking as t — )■ 0~ subsets will eventually help to better 
understand the possible nature of other plausible blow-up singularities of the NSEs. As a 
common feature of our blow-up analysis, we will see that a complete/closed set of vector 
solenoidal Hermite polynomialqj in M'^ can play an important part. 

Our main goal here is as follows: using techniques of blow-up scaling and matched 
asymptotic expansions from react ion- diffusion theory, to show that a Petrovskii 's-like cri- 
terion of boundary regularity for the heat equation (1934) does occur for the NSEs (II. ip . 

1.3. Petrovskii's criterion of 1934 for the heat equation: a first discussion. For 

the heat equation, 

(1.6) ut = Au in Qo, 

the regularity problem of the characteristic vertex (0,0) was optimally solved for dimen- 
sions A^ = 1 and 2 by Ivan Georgievich Petrovskii in 1934-35 [^ 16 Ip . who introduced 



These are the same vector polynomials that occur in the study of multiple zero formation for the 
Stokes equations and NSEs, [HI § 3]. 

Compare with the dates of Leray's pioneering study, "1933-34" ; almost a perfect coincidence! 
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his famous Petrovskii 's regularity criterion; see [T7] and [20] for a full history and further 
developments in general parabohc theory. This is the so-called "2v^log log-criterion" (see 
(12. 8p and (12. 9p below), which we are going to achieve, at least formally, for the 3D NSEs. 

The following issues naturally occur: 

(i) On one hand, Petrovskii's-like criterion can be expected, since (II. ip . similar to (II. 6p . 
is indeed a (vector) parabolic second-order equation; 

(ii) On the other hand, (II. ip is a nonlocal parabolic PDE for solenoidal vector fields, 
and it is not straightforward from the beginning that this cannot affect the regularity 
standing; and 

(iii) Finally and most essentially, (II. ip contains both linear (the second-order Laplacian) 
and nonlinear (convective) operators, so that the regularity of the vertex (0,0) inevitably 
will depend on both, which makes the analysis more difficult. 

Note that both the issues (i) and (ii) apply to the linear Stokes problem without the 
quadratic convection (II. 4p . for which our regularity results turn out to be new as well. 

1.4. Layout of the paper. In Section HI we perform first a blow-up scaling near the 
characteristic vertex (0,0). In Sections E] and IH the main goal is to show how the convec- 
tion term in the NSEs (II. ip can affect the regularity conditions by deriving sharp formal 
asymptotics of solutions near the characteristic point. A necessary and already well- 
existing spectral theory involving a complete set of vector solenoidal Hermite polynomials 
as eigenf unctions of the linear Hermite operator is described in Appendix A. 

For the sake of our regularity study, we apply a method of a matched asymptotic 
(blow-up) expansion, where a Boundary Layer behaviour close to the lateral boundary 
OQq (Section [3]) is matched, as t — )■ 0~, with a ''centre subspace behaviour" in an Inner 
Region, developed in Section |H This leads to a perturbed 3D nonlinear dynamical system 
for the first Fourier-like coefficients in the eigenfunction expansions of the vector field 
u(x, t) via standard solenoidal Hermite polynomials. This approach falls into the scope 
of typical ideas of asymptotic PDE theory, which got a full mathematical justification for 
many problems of interest. We refer to a most general asymptotic analysis performed in 
and also to a number of complicated blow-up asymptotics in reaction-diffusion theory 

|. According to the classification in [31], our matched blow-up approach corresponds to 
perturbed three-dimensional dynamical systems, i.e., to a rather not-that-advanced case 
being however a constructive one that has given a number of new asymptotic/regularity 
results. We propose a final, more general discussion of various (boundary and interior) 
blow-up singularities for the NSEs in Section 

In addition, for more clear expressing our regularity techniques and their applicability 
in general PDE theory, we develop in Appendix B (the C one contains the corresponding 
Hermite spectral analysis) at the paper end, as a natural extension, a similar regularity 
analysis of the well-posed Burnett equations in Qo with zero Dirichlet conditions 

fUi + (u ■ V)u = — Vp — A^u, divu = in Qq, 
u = Vu ■ n = on dQ^, u(z, — 1) = uo(a:), 



where n denotes the unit inward normal vector to dQor\{t}. Here we have the bi-harmonic 
diffusion operator — A^u on the right-hand side of the u-equation. It turns out that our 
general scheme of the boundary regularity analysis can be applied; however, with harder 
asymptotics and more formal nature of the final difficult estimates. 

Concerning other problems and techniques of modern regularity theory, we refer to 
monographs [271 IMl EH [51], |5l] and [35], [l6]-[53] as an update guide to elliptic regularity 
theory including higher-order equations, as well as to references/results in |28l [381 SD US 
69] and [T71 [20] for linear and semilinear parabolic PDEs. 



2. First blow-up scaling: Sturm's backward scaling variable, paraboloid 

GEOMETRY, AND THE CaUCHY PROBLEM 

2.1. First blow-up scaling: exponentially small perturbations of a rescaled par- 
abolic flow. We perform blow-up scaling in (11.11) for the regularity analysis 

(2.1) u(x,t)=v(2/,T), ?/=^, r = -ln(-t):(-l,0)^M+, 

where y is, indeed, Sturm's backward rescaled variable introduced him in 1836 [65j in the 
study of zero sets of solutions of linear parabolic equations such as (II. 6p for A^ = 1. 
Thus, scaling (12. ip yields the following exponentially perturbed rescaled equatioiQ: 

Vt- + e~2 (v ■ V)v = — e~2 Vj9 + B*v, divv = in Qo, 
where v = [v-^, f^, v^]'^ and B* = A — | y ■ V 
is Hermite's classic symmetric (self-adjoint) operator, [H p. 48]. 

2.2. Backward paraboloid geometry and a slow growing factor (/^(r). According 
to Petrovskii [501 EI]) the backward paraboloid dQo will be defined as follows: it is a 
perturbation of the standard fundamental backward one. 



(2.3) S{t) = dQo n {t} : go(^) = ^JEll (^^A = H)'^ ^(r), 

where Oj G (0, 1] are normalized constants. We can treat more general cont'ex paraboloids, 
but, for simplicity, will restrict to the basic ones as in (12. 3p . which are also of a challenge. 
For difficult estimates to follow, we will use radially symmetric paraboloids with Oj = 1: 

(2.4) S{t) = dQo n {t} : qo{x) = \x\ = (-t)l </.(r). 
In (12. 3 p and (12. 4p . (/^(r) > is a slow growing function satisfying 

(2.5) (/?(r) -^ +00, if'ir) > 0, <^'{t) ^0, and ^ -> as r -^ +oo. 

Moreover, as a sharper characterization of the above class of slow growing functions, we 
use the following criterion: 

(2.6) (^)'^oo as r^+oo. 



Here, (|2.2[) is not the rescaled one (15.31) written in Leray's variable (|5.2p . a difference to be discussed. 
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This is a typical condition in blow-up analysis distinguishing classes of exponential (the 
limit is 0), power-like (a constant 7^ 0), and slow-growing functions. See [HSl PP- 390-400], 
where in Lemma 1 on p. 400, extra properties of slow-growing functions f l2.6p are proved. 
For instance, one can derive the following comparison of such f{T) with any power: 

(2.7) for any a > 0, (/?(r) < r" and ip'ir) < r""^ for r > 1. 

Such estimates are useful in evaluating perturbation terms in the rescaled equations. 

In Petrovskii's criterion for the heat equation (II. 6p . for any A^ > 1, the "almost optimal" 
function, satisfying (12. 5p . (12. 6 p and delivering a regular vertex (0,0), is 



(2.8) V^*{t) = 2y/hiT as r -> +cx). 



Replacing this fundamental constant "2" by 2-l-e, for an arbitrarily small constant e > 0, 
makes (0,0) irregular for the heat equation (II. 6p . In the general case of arbitrary (p{t), 
Petrovskii's criterion, for our A^ = 3, for the radially symmetric paraboloid (12. 4p readq3 



(2.9) 



/^3(^)e-^^W/4dr 



00. 



These dependencies will be compared with those obtained for the NSEs (11.11) . 

Thus, the monotone positive function v'(r) in (12. 3 p is assumed to determine a sharp 
behaviour of the boundary of Qq near the vertex (0,0^) to guarantee its regularity. It 
follows that the rescaled equation (12. 2 p is set in an expanding rescaled domain 



(2.10) S{r) = dQo n {r} : qo{y) = ^JEi=l ^iVf = ^{r) ^ +00 as r ^ +00. 

By n we denote the inward unit normal to S{t). In the limit r = -fcxo, we arrive at the 
equation (12. 2p in the whole space M.^, which requires some spectral theory (Appendix A). 

2.3. Towards the Cauchy problem. In Inner Region (see Section [3] for details), de- 
scribed by compact subsets in the variable y in (12.11) . we deal with the original rescaled 
problem (12. 2p in the unboundently expanding domains (I2.10p . As usual and customary 
in potential and general PDE theory (see e.g., Vladimirov [68j), it is convenient to con- 
sider the NSEs in whole space M'^ x [—1,0). Note that, in the study of the NSEs in 
non-cylindrical domains, Fujita-Sauer [T5l[T6] also extended the problem to M^ by intro- 
ducing a strong absorption term —nu in the complementary domain on the right-hand 
side in (II. ip and passing to the limit n — )■ +00. Then, in view of the control of the 
total kinetic L^-energy, this "regularized" solution u = u'-"-' — > as n — > +00 outside the 
given non-cylindrical domain, and hence, in the limit, the zero Dirichlet conditions on the 
boundary are restored. 

Thus, we extend v(?y,r) by beyond the boundary, i.e., set 

,2.11) .fo..,^vte..,^,.M-,„,.))^{^<'";' I: °*7y)>l(rf'' 



^For Qq C K^ X [—1,0), the multiplier ip"^ in this integral criterion is replaced by ip^ . 



where H{-) is the Heaviside function. Then, since v = on the lateral boundary S{t) = 
{Qo{y) = '^{t)}^ one can check that, in the sense of distributions (see e.g., [68, § 6.5]), 

(2.12) ^r = ^rH, Vv = Vv//, Av = Avi7+g(5^(^), 

where a single-layer potential with the density /W = |^ acts as follows: for any G C^, 

(2.13) {f^Ssir)^(P) = Is(r) /^'^d^- 

In order to avoid a pressure trace on S{r), we perform a continuous extension oi p{y, r) 
by solving the Laplace equation in the outer domain: 



(2.14) Ap = in R^ \ (Qo(t-) n {r}), p = p on S{t). 

This outer Dirichlet problem is known to admit a unique solution p{y, r) vanishing at 
infinity, [681 § 28], so that we are given the continuous pressure piy^r) defined in the 
whole M^ X R I . It then follows that 



Vp = Vp in Qq{t) n {r} (since the jump is zero: [p\s(t) = 0)' ^^^^ 
(2.15) 

Vp is div-free in R^ as a distribution : (Vp, </>) = 0, V0 G C^, div = 0. 
Thus, {v,p} satisfies the following Cauchy problem in R^ x R_,_: 

(2.16) v^ + e-i(v- V)v = -e-5 Vp + B*v-g%(^), divv = 0. 

Hence, we obtain a single perturbation term on the right-hand side expressed in terms of 
a simple layer potential with the prescribed density on the surface (I2.10p . changing with 
the time r. Clearly, various linear and nonlinear "interactions" of all these and other 
operators in (I2.16P will define regularity of the vertex. 

2.4. The Cauchy problem in Leray's nonlocal setting. Using Leray's nonlocal for- 
mulation of NSEs [45i p. 32], we next apply to the Cauchy problem (I2.16P the operator 

P = /-VA-i(V-/) (||P|| = 1) 

being the Leray-Hopf projector of (L^(R^))^ onto the subspace {w G {L'^Y '■ div w = 0} 
of solenoidal vector fields. Let us note another representation of the projector P therein: 

Pw = \vi — i?icr, V2 — -R2C", V'i — R-ifyY , where a = RiWi + R2W2 + RsW^, 

and Rj are the Riesz transforms, with symbols fe. Using the fundamental solution of A 
in R^, A^ > 3, and denoting by a^ the surface area of the unit ball Bi C R^, 

(2.17) bNiy) = - (jv-Viv M^' ^^^^^ ^^ = mm ^^3 = 47r), 
the projection of the convective term reads: 

-P(v . V)v = -(v • V)v + C3 / r^ tr(Vv(^, r)f dz, 

(2.18) K3 

where tT{V^r{z,T)f = ^i^,) vl^v^ and C^ = -L > q (C3 = ^). 

This is a more convenient form for some estimates. 
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Using the projection P eliminates the pressure term V^ in f l2.16p . and we obtain the 
Cauchy problem for the following perturbed nonlocal parabolic equation for v: 



(2.19) 



B*v-e-^P(v- V)v-g(5s(^) in R^ x 



Since, by construction, the last term is solenoidal, we have omitted the projection P 
therein. We recall again that local existence and uniqueness of a classic solution v(?y, r) 
of fl2.19p are guaranteed by known local regularity properties of the NSEs. Moreover, for 
any sufficiently small data vq, solutions of f l2.19p are well defined for all r G M+, i.e., up 
to the boundary blow-up moment t = 0~ (t = +cxd). For other solutions, in general, we 
assume that v(?/,r) are well defined and do not blow-up at a finite r > 0, so we need to 
study their behaviour as r — ;■ +00. 

It then follows that Wiener's regularity of the vertex (0, 0) is equivalent to the following: 



(2.20) 



is globally asymptotically stable for (12.191) . i.e., 

any such solution of (I2.19P satisfies v(?/, r) — )■ as r — )■ +00 uniformly in y. 



2.5. A full pressure representation in R^ x M+. As usual [l5l p. 30], once the vector 
field V has been obtained from (I2.19p . the pressure is then given by the corresponding 
Poisson equation. Since it is slightly more technical to get it from the rescaled equation 
(I2.19P containing extra operators, we consider first the original (non-rescaled) problem for 
u = uif (■), where the Heaviside function H is concentrated on Qo{t) fl {t}, and construct 
a harmonic extension p as in ( 12.14p for S{t). This equation, which will be also in use, is 

(2.21) Ui + (u-V)u = -Vp + Au-|i55(t), 

where we keep the same notation and arguments as in (12.190 . Then, taking div [45, p. 30], 
we obtain the following equation with two extra densities of special potentials: 

(2.22) -A,p = tr (V,u)2 + div, {^ Ssn)) - [V.p]s(t) Ssn)- 

Observe that the jump of the pressure gradient [Va;p]5(j) enters the second density, that 
makes this elliptic problem more nonlocal and hence more difficult. 

After scaling in (12. ip . i.e., setting x = yy/—i when approaching the vertex (0,0), we 
have from (I2.22p . taking into account that 6s(t) = i^t) ^s(t)^ 



(2.23) -Ayp = tr (V,v)2 + e- div, (£ 6^..) - e"^ [V,p]^. . S, 



m 



Therefore, this nonlocal problem for p is reduced to a Fredholm linear integral equation 
(with a positive kernel) of the second kind, 

(2.24) p = (-A,)-i [tr (V,v)2 + e^ div, {^ 5^^^^) - e"^ [V,p]^(,) 6^^^^] , 

where (— A,)^^ is defined by the convolution with the fundamental solution (I2.17p . A = 3. 
It follows that the behaviour of p as r — )■ +00 is dependent on the trace of its gradient 
on the expanded boundary S{t). Fortunately, as r — )■ +00, the jump of the gradient 
[Vyp]^/^-) on the right-hand side of (12.240 has an exponentially small infiuence on p, so 
that (I2.24p /( l2.23p are "almost" standard integral/elliptic equations with good positive 



kernels/Laplacian. However, we will need to justify that, nevertheless, the corresponding 
densities of these single and "double-layer-type" potentials do not get exponentially large, 
thus undermining their exponentially small factors in front of them. This can be done a 
posteriori, when the independent rescaled v-problem (12.191) has been solved. 

On the other hand, introducing in (12.231) the following integral operatoo 

(2.25) M(r)p = -Ayp + e"^ ['^yP]sir)^sir), 
the pressure is given by 

(2.26) p = M-i(r) (tr (V.v)^ + e^ div, (|j S^^^^)) . 

Indeed, this is the equivalent pressure representation, since M^^(r) is well defined by 
the local well-posedness of the NSEs in bounded domains with smooth non-characteristic 
boundaries. 

3. Boundary layer expansion close to dQo 

3.1. Two region expansion. As we have mentioned, by the divergence of fir) in (12. 5p . 
sharp asymptotics of solutions close to the vertex (0, 0) will essentially depend on the 
spectral properties of the linear operator B* in the whole space R^ (see Appendix A), as 
well as on the nonlinear projected convective term. This "interaction" between linear and 
nonlinear operators in the NSEs, together with the paraboloid shape, are key for us. 

Studying asymptotics of solutions of the rescaled problem (12. 2p . as rather often occurs 
in difficult blow-up expansions in nonlinear PDE theory, this singularity problem is solved 
by matching of expansions of solutions in two regions: 

(i) In an Inner Region, which is situated around the origin y = 0, and 

(ii) In a Boundary Region close to the boundary (I2.10p . where a Boundary Layer occurs. 

In other words, we show that generic behaviour of solutions of the NSEs in shrinking 
neighbourhoods of the paraboloid vertex (0,0) is not of any self-similar form, and hence 
gets more complicated and demands novel non-group-similarity techniques to detect. 

Actually, such a two-region structure (i)-(ii) above, with the asymptotics specified 
below, defines the class of generic solutions under consideration. We begin with a simpler 
analysis in the Boundary Region (ii). 

3.2. Boundary layer (BL) variables and perturbed equation. Sufficiently close to 
the lateral boundary of Qq, it is natural to introduce the next rescaled variables 

(3.1) ^=^ ^^^ v(|/,r) = w(z,r). 

This makes the corresponding rescaled paraboloid (I2.10p fixed: 



(3.2) S: J^=,a.zf = l. 



^3 



For T ^ 1, this is just an asymptotically small perturbation of the Laplacian; though proving that p 
on S{t) does not grow exponentially fast is a part of the problem. 
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The rescaled vector field w now solves a perturbed equation: 

(3.3) w. = ^A,w-i^-V,w + ^^-V,w-ie-iP(w-V,)w-(^|^,5^. 

Let us introduce the BL- variables: fixing a point ^q G S* on the boundary (13.21) . we set 

(3.4) i = ip'^{T){zQ- z) = ip{ipzQ-y), ip'^{T)dT = ds, w(2;, r) = p(s)g(^, s), 

where p{s) = [p^(s), p^(s),p^(s)]^ G M.^ for s ^ 1 is an unknown scaling slow vary- 
ing/decaying (in the same natural sense, associated with (12. 6p ) time- factor depending on 
the function (p{t). Cf. e.g., ~ -rp: as a clue. As usual, this p-scaling is chosen to get 
uniformly bounded rescaled solutions, so, we naturally assume that, for each component, 

(3.5) supg g^{^, s) = l for all s > 1, j = 1, 2, 3. 

On substitution into the PDE in (13. 3p . we obtain the following perturbation of a linear 
uniformly parabolic equation: 

_£ e-i P(g ■ \/^)g - ^ 1^ Ss^ir), where Ag = A^g + i ^o ■ V^g, 
and S^{t) is the boundary (13. 2 p expressed in terms of the BL- variable ^ in (13. 4p . so that 

^•55 (r) = ^3^ ^S- 

As usual in boundary layer theory, the BL-scaling (13. 4p means that we are looking for 
a generic pattern of the behaviour described by the perturbed equation (13. 6 p on compact 
subsets, shrinking (focusing) to a fixed Zq on the lateral boundary, 

(3-7) \^\ = o{^)^0 =^ \z - zo\ = o{^^) ^ as r ^ +oo. 

Thus, in (13. 6p . we arrive at a linear uniformly parabolic equation perturbed by a number 
of linear and nonlinear terms, which, under given and other special hypothesis to be 
specified, are asymptotically small. Indeed, on the space-time compact subsets (13. 7p . the 
second term on the right-hand side of (13.60 becomes asymptotically small, while all the 
other linear ones even smaller in view of the slow growth/decay assumptions such as 
(12. 6p for V9(r) and p{s). In particular, the rescaled nonlinear convective term in (13. 6p 
is asymptotically small on bounded rescaled vector fields g in view of an exponentially 
decaying factor and by the hypotheses (12.50 . 

However, the last term in (13. 6p . given by a density of a simple layer potential, requires a 
special treatment. Indeed, this density depends upon a still unknown gradient of w = pg 
on the boundary, which is under scrutiny in the present BL-analysis. However, the nature 
of our BL-scaling assumes that we deal with uniformly bounded rescaled function g, on 
which the last term is asymptotically small and is of order ~ -4:p- — )■ as r — )■ -|-oo. 

The BL-representation (13.40 . by using the rescaling and (13. 5p . naturally leads to the 
following asymptotic behaviour at infinity: 

(3.8) lims-^+oo5'^(^, s) -^ 1 as ^ -^ oo, 
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where all the derivatives also vanish by the standard interior parabolic regularity. Actu- 
ally, the nature of the BL-scaling (13 ■4p near the point Zq E S implies that, asymptotically, 
the limit problem becomes one- dimensional, depending on the space variable 

(3.9) ^ = e-n, 

where n is the unit inward normal to the smooth boundary S in fl3.2p . Therefore, the 
limit ^ — )■ oo in fl3.8p should be also understood in the sense of the single variable fl3.9p . 
This essentially simplifies the BL-structure to appear. 

Moreover, since according to the BL-scaling (13. 7p . as s — t- +oo, the rescaled solution 
becomes constant (see (13. 8p ) a.e. and hence solenoidal, we do not need to require the limit 
BL-profile to be solenoidal as well. Moreover, we will see that, in the actual boundary 
layer, the BL-asymptotic is "almost" solenoidal, up to an exponentially small perturba- 
tion. 

3.3. Passing to the limit: generic solutions. Thus, we arrive at the problem of 

passing to the limit as s — t- +oo in the problem (13. 6p . (13. 8p . Since, by the definition in 
(13. 4p . the rescaled orbit {g(s), s > 0} is uniformly bounded, by classic parabolic interior 
theory [9l [TOl [11], one can pass to the limit in (13. 6p along a subsequence {s^} — )■ +cxd. 
Namely, we have that, uniformly on compact subsets defined in (13.70 . as A; — )■ oo, 

(3.10) g(sjt + s) — )■ h(s), where h^ = Ah, h = at r] = 0, h^rf=+oo = '^■ 
Consider this one-dimensional limit (at s = +oo) equation obtained from (13. 6p : 

(3.11) h, = Ah = h^^ + |h^ in M+ x M+, h(0,s)=0, h^{+oo,s) = l. 

It is a linear parabolic PDE in the unbounded domain IR+, governed by the operator A 
admitting a standard symmetric representation in a weighted space. Namely, we have: 

Proposition 3.1. (i) (13. lip is a gradient system in a weighted {LP'Y-space, and 

(ii) for hounded orbits, the uj-limit set VLq of (13.110 consists of a unique stationary profile 

(3.12) gl[^) = I - e~^l\ J = 1,2,3, 

and VIq is uniformly stable in the Lyapunov sense in a weighted {L'^)^ -space. 
Proof, (i) As a second-order equation, (13.110 can be written in the symmetric form 

(3.13) e'^/^K = (e^/'h,),, 

and hence admits multiplication by h^ in (L^)^ that yields a monotone Lyapunov function: 

(3.14) me^^\hr)^ = -je^/\hif<0. 

Note that, regardless that the weight e'''^ is exponentially growing as 77 — )■ +00, on the 
limit profile (13.120 . all the functionals in (13.140 are well defined. In other words, the 
problem (13.60 is a perturbed gradient system, that makes much easier to pass to the hmit 
s — !■ +00 by using power tools of gradient system theory; see e.g.. Hale [29] . 
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(ii) For a given bounded orbit {h{s)}, denote h^{s) = Qq + w^^s), so that w(s) solves 
the same equation fl3.13p . Multiplying by w(s) in {L'^Y yields 

(3.15) l£ J e'^/^w^f dr] = - J e'^/^w^^f dr] < 

for any nontrivial solution, whence the uniform stability (contractivity) property. D 

Finally, we state the main stabilization result in the boundary layer, which also estab- 
lishes the actual class of generic solutions we are dealing with. 

Proposition 3.2. Under specified above assumptions and hypotheses, there exists a class 
of solutions of the perturbed equation f l3.6p . for which, in a weighted {L'^)^-space and 
uniformly on compact subsets, 

(3.16) 9'{C,s)^gl{0 as s^+oc (j = 1,2,3). 

Proof, (i) Under given hypotheses, the uniform stability result in (ii) of Proposition 13.11 
implies [231 Ch. 1] that the w-limit set of the asymptotically perturbed equation (13. 6p is 
contained in that for the limit one (13. lip , which, under the given hypotheses, consists of 
the unique profile (I3.12p . D 

3.4. BL-behaviour is "almost" divergence free. For the future convenience, we state 
again the asymptotic BL-behaviour: in the rescaled sense, by (13.160 . (I3.12p . and (14. 8p . 



(3.17) \r{y,T) = Co{T)go{y,T) + ..., where 5-^(7/, r) = 1 - e" 



^dist{y,dQo{r)} 



It is important that, by (I3.17p . the first term is "a.e." an exponentially small perturbation 
of a divergence-free flow. Indeed, differentiating (13.170 at any point staying away from 
the boundary by an arbitrarily small fixed dist{-} = 60 > 0, we have 

(3.18) divj,v(y,r) = 0(co(r)(^(r)e-^'^«) ->0 as r -> +cx), 

provided that Co(r) is not exponentially large (this does not happen, as we will show). 
In other words, not that surprisingly, the BL-expansion well keeps solenoidal features of 
originally divergence-free solutions v and, as customary, just makes an asymptotically 
(i.e., exponentially as in (I3.18P ) small perturbation of the div. Thus, a somehow essential 
violation of the solenoidal property can happen only in an asymptotically small 0(-7:pr)- 
neighbourhood of the boundary, which is negligible and plays no role for r ^ 1. 

4. Inner Region expansion: towards ODEs regularity criterion 

4.1. A standard semigroup approach leads to a more complicated problem. Let 

us first perform necessary manipulations using a standard semigroup approach. Applying 
to (I2.2ip the projector P yields 

(4.1) Ut = Au-fH<55(t)-P(u-V)u. 
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Therefore, using the fundamental solution b{x, t) of the heat operator Df — A with the 
rescaled kernel (the Gaussian) as in flA.Sp gives the following convolution representation 
of the solution of (14. ip : 

(4.2) u(t) = b{t) *uo- J b{t-s)* g(s) 5sis) ds-Jb{t-s)*F (u(s) ■ V)u(5) ds. 



In particular, taking div, we see that the second term on the right-hand side, containing 
a surface integral, is div-free, since u is. Finally, sharply estimating the normal derivative 
therein and in the third term from the core of BL-theory, the asymptotics (I3.17p . one can 
study the asymptotic behaviour of solutions, after using the necessary scaling (12. ip . 

However, it turns out that this standard integral semigroup approach leads to a more 
complicated analysis, than a differential one we will perform by using known spectral 
properties of the rescaled operator B* involved in (I2.19p . Nevertheless, it is worth men- 
tioning that such an approach can be translated to the integral equation (14. 2p . with the 
clear advantage of a more reliable rigorous justification by using obviously smoother prop- 
erties of solutions and, as a result, their better uniform estimates in stronger metrics. It 
should be noted that some principle difficulties cannot be avoided in such a way, and, 
overall, technical technical questions become more dominant for (14. 2p . 

4.2. Eigenfunction expansion: derivation of a 3D dynamical system. Thus, in 
the Inner Region, we deal with the original rescaled Cauchy problem (I2.19p . Since, by 
construction, the extended solution orbit (12. lip is uniformly bounded in (L^*(M^))^, we 
can use the converging in the mean (and uniformly on compact subsets in y) eigenfunction 
expansion via the solenoidal Hermite polynomials as in flA.17p : 

(4.3) v(2/,r) = E(;3)C/3(r)v;(2/), 

where we use the convention (lA.lSp . Substituting (14. 3 p into (I2.16P and using the or- 
thonormality of these polynomials yield the following dynamical system for the expansion 
coefficients: for all |/3| > 0, 

(4.4) tp = Xpcp - (^(^^M' v^) - e-5 (P(v ■ V,)v//, v^,), 

where Xp = —-^ ^i-re the real eigenvalues as in (lA.2p . 

Recalling that eigenvalues in (14. 4p satisfy Xp < — | for all \j3\ > 1, it follows that we 
need to concentrate on the "maximal" first Fourier generic pattern associated with the 
first constant Hermite polynomial Vq in (lA.14p . 

(4.5) k=\(3\ = 0: Ao = and v*(y) = e = [1, 1, 1]^. 
The normalized eigenfunction of the L^-adjoint operator B is then 

(4.6) My) = F{y)e, 

where F{y) is the Gaussian in ( IA.6p . Actually, as follows from (I3.16p . this corresponds to 
a naturally understood "centre subspace behaviour" for the operator B* in (I2.16p : 

(4.7) ir{y, r) = Co(r)e + w^{y, t), where w^ G Span{v;, |/3| > 1}, 
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and w^{y, r) is then negligible relative to Cq{t) for r ^ 1. 

Fortunately, we actually do not need such a literal using of those "centre subspace 
issues" for a difficult non-autonomous equation like f l2.19p . since the asymptotics of so- 
lutions (14. 7p is directly dictated by Proposition 13.21 In its turn, this is an equivalent 
characterization of our class of generic patterns, and, in particular, the following holds: 

Proposition 4.1. Under the given above assumptions and hypotheses, for the prescribed 
class of generic solutions defined in Proposition \3.2. (14.71) holds with w-^^{t) = o(|c^(t)|) 
as T ^ +00, and then the matching with the boundary layer behaviour in (13.41) requires 

(4.8) ^^1 as T^+oo =^ f>i{s) = ai{T){l + o{l)), j = 1,2,3. 

Proof. This follows from the construction of the boundary layer by comparing the solution 
representations in (13. 4p . (I3.16p . and (14. 7p . D 

Thus, the equation for Cq{t), with Aq = 0, takes the form: 

(4.9) Co = -(£<55M' vo) - e-i (P(v ■ V,)vi/, vq), 



where the first adjoint eigenfunction vq is as in (14. 6p . 

We now need to return to BL-theory in Section [3] establishing the boundary behaviour 
(13. 4p for T ^ 1, which has the form (I3.17p . Then the convergence (I3.17p . which by a 
standard parabohc regularity is also true for the spatial derivatives, yields, in the natural 
rescaled sense, 

(4.10) £ = -| Co(T)(y9(r)^^ + ... , where y ■ n{y) < (by convexity). 

Therefore, the first (linear) term in (14. 9 p asymptotically reads 

(All) _/9v^ \ ^ _co^(iL r s-"Wp-kiV4^ I 

y^-'-'-) \an"5(T)'^0/ 2{47r)3/2 J \s\ ^ aS -\- ... . 

S{r) 

Due to the normalization in (12. 3p . we have that Oj < 1, so that the last term can be 
estimates above as: for any component j = 1, 2, 3, 

(4.12) ^^ / ^e-NV4d5<-7iV^^(r)e-^^W/^ + ..., 71 = const. > 0. 

S{t) 

Recall that the extra factor tp"^ is obtained via integration over a closed surface in M'^. 
For the simple radial paraboloid shape in (I2.10p and (13. 2p . with a^ = 1, i.e., by (12. 4p . 

(4.13) Sir)=dQon{r}: qoiy) = \y\ = ifir), and S : \z\ = 1, 

(I4.12P presents a sharp asymptotics behaviour rather than an estimate: for j = 1, 2, 3, 

(4.14) 0^ J ^ e-l^lV4 d. = -Co7i¥^=^(r) e'-^M/^ + ... . 

S{t) 

Naturally, it is possible to derive most sharp asymptotics for this radial case. 
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Finally, we need to estimate the last nonlinear quadratic term in f l4.9p : by f l3.17p . 

(4.15) e-i (P(v ■ V,)v//, vo) = e-5 (cq ■ e)co ^^ / P[(go ■ ^y)goH] e'l^l'/^ dy + ... , 

where the quadratic term (cq ■ e)co denotes the vector (cq + Cq + Cq)co G R^. Fortunately, 
since this nonlinear term enjoys having a fast decaying exponential factor e~2, we do 
not need its better sharp estimates. We just need to show that, on the generic solutions 
obeying the BL-behaviour ( 13.1 7p . this term always remains exponentially small, so does 
not play any role for the regularity conclusions. 

We restrict to the radial case, though the asymptotic smallness of the convection can 
be similarly shown for more general convex paraboloids. Thus, using the representation 
of P(v ■ V)v given in fl2.18p and using the change as in (13. ip . i.e., setting y = '^{t)z, we 
have from (I3.17p . for r ::^ 1, 



^2 



e 2 (v ■ V)v = e 2 (cq ■ e)co V5[(l — e 2 '^^^e 2'^^(n-e) 
-C,^' I ^ e-^^'^c ^^^^^.^(n . e,)(n ■ e,) dC] + ... = Ji + J2, 

where d^ (and (i^ in the integral) denotes the distance: d^ = dist {z, S}. For such a rough 
estimate from above, one can omit the projector P, using the fact that ||P|| = 1. 

It is not difficult to see that, the Ji-term in (I4.16P is leading to the following integral: 

2 2 2 

(4.17) Ji ~ e"^ (co ■ e)co(/5^ / e"^ l^^l (l - e"^'^")e"^"'"(n ■ e) dz, 

{kl<i} 

where, in the radial geometry, we may put d^ = 1 — \z\. Reducing the integral in (14.171) 
to a standard ID one, it can be estimated above as follows: for some constant 72 > 0, 

2, X 

(4.18) |-^i(i")| < 72©"^ |(co ■ e)co| v5^(T)e"^^^ as r -)■ +00. 

Consider the last term J2 in (I4.16p . The corresponding upper estimate is: for 734 > 0, 

l^2(r)|<73e-^|(co-e)co|v^«| / e"^ N^ ^ e'-^'^c ^^^^.(.) dCd.| 

(4.19) il-l<i} 

< 74e 2 |(co ■ e)Co\ (p^{T)e 4 . 

We do not guarantee that the 99^ multiplier in the final estimate in (I4.19P is any sharp 
(as well as y?^ in (I4.18P ). but it is sufficient for showing the convection neglect near the 
vertex. Indeed, any such very rough estimates (or omitting the projector P as we did 
above) cannot undermine the principal fact: extra multipliers containing any power of 
(/?(r), being a slow growing function, do not practically affect the exponentially decaying 
factor e"'^/2 as r -;■ +00 in (I4.15P and (I4.16p . Comparing with (I4.17p . (I4.18P yields that 
(I4.19P supplies us with the leading term as r — )• +00. 

Thus, bearing in mind all above assumptions and estimates for generic patterns, we 
obtain the following asymptotically approximate dynamical system for the first expansion 
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coefficients Cq{t): for r ^ 1, 



(4.20) Co ~ -7i Co^^t) e-^'M/4 + 74 e"^ (cq ■ e)co ip%T) e-'^'M/^ + ... 



where we now omit all higher-order terms appeared via the above hypotheses. The sign 
"~" in fl4.20p means that, in the presentation of the influence of the nonlinear convec- 
tion term, we used the estimate fl4.19p . rather than a sharp asymptotics. However, this 
estimate suffices to declare that the convection term is negligible in the regularity analysis. 

4.3. 3D regularity criterion. Thus, according to (12.201) . the following conclusion holds: 

Theorem 4.1. Under the assumed above hypotheses and conditions, the vertex (0,0) is 
regular for the NSEs problem (II. ip . (I1.2p in the class of generic solutions, iff 



(4.21) the origin is globally asymptotically stable for the 3D dynamical system {4-20) 



i.e., any its solution satisfies: 

(4.22) Co(t) ^0 as r -> +cx). 

4.4. Two regularity conclusions. We begin with a simpler linear one. 

1. Linear Stokes problem . As a ffist consequence, we con&m that Petrovskii's cri- 
terion (12. 9 p remains valid in the linear case. Recall that here, our analysis do not include 
more difficult "nonlinear" estimates used in (I4.15p . 

Theorem 4.2. For the linear Stokes problem (II. 4p . under our hypotheses on generic 
solutions, the regularity criterion of the vertex (0,0) is Petrovskii's one (12. 9p . 

Proof. Introducing the new time, 

(4.23) (/?=^(r)e-^'W/^dr = ds ^ s = J ip^iQe-^''^^'^/^ dC, 



from (I4.20p (with no quadratic term), we obtain a linear diagonal autonomous system 

(4.24) A Co = -71 Co =^ Co{s) = Co(0) e"^^^ -> 

if and only if s — > -|-oo as r — )• -|-oo, and the divergence in (12. 9p follows. D 

2. Navier-Stokes equations . Performing the change (I4.23P for the full dynamical 
system (I4.20p yields 

(4.25) £^0 = -7iCo + 74e"5 (cq ■ e)co (/^^(r). 

An elementary balancing of the linear and nonlinear term on the right-hand side shows 
that the nonlinear one can be efficiently involved into the regular asymptotics if it has at 
least an exponential growth 

(4.26) |co(r)|»ei;^»l for r » 1, 

since, by assumptions, (/?(t) is a slow growing function. Of course, (I4.26P is impossible, 
since by the regularity assumption Co(t) -^ 0. 
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Note that, if v(r) gets vanishing as r — t- +oo, then the same is true for the pressure 
via f|2.26p (then f l2.25p is a perturbation of the Laplacian), i.e., the pressure influence in 
the full equation (with no convection) is truly exponentially negligible. 

Thus, this is our final conclusion for the NSEs: under given hypothesis and for generic 
solutionqj, the nonlinear convection term cannot affect the regularity of the vertex (0,0), 
so that Petrovskii 's criterion ( 12. 9^ remains true for the Navier-Stokes equations (II. ip . 



5. Final discussion: two key blow-up problems for the 3D NSEs 

It is now worth and natural to look how the characteristic boundary regularity analysis 
stands and fits in the wide area around the Millennium Prize Problem for the Clay Insti- 
tute (the MPPCI) on global existence or nonexistence of bounded smooth solutions for the 
NSEs (II. ip : see Fefferman |14j. Both deal with settings presented in a similar fashion: 

Boundary point Regularity: sharp "blow-up" asymptotic expansions of pos- 



sible solutions of (II. ip near a characteristic boundary point of Qq. Then, mathematically 
speaking, after blow-up scaling (12. ip . we arrive at an exponentially small perturbation of a 
solenoidal heat equation in (12. 2p . which is convenient to write down in its full presentation 



(5.1) 



v^ = Av — \y- Vv — 


T 

e 2 


P(vV)v in Qo- 



We observe here a crucial moment: in the vertex regularity study, the key nonlinear 
convection term gets an extra fast exponentially decaying factor e~2 (boxed in (15. ip ). 
that, indeed, essentially simplifies the matched asymptotic analysis. As we have seen, 
our goal was then simply to show that this nonlinear exponentially small perturbation 
does not affect the regularity of the vertex (0, 0) at all, so it becomes pure parabolic 
(Petrovskii's) one, i.e., governed as r — t- +oo by the rescaled solenoidal heat equation 

v^ = Av - i y ■ Vv in Qo- 

(IR) Interior point Regularity (the MPPCI): asymptotic expansion of possible 
solutions in any interior point to check whether finite time blow-up in L°° is possible or 
not. Then, a full Leray's self-similar blow-up scaling must be performeci 

(5.2) u(x,t) = ^v(2/,r), 1/=^, r = -ln(r-t). 



For the full Navier-Stokes models, in view of an essential difficulties to justify nonlinear convection 
estimates for calculating (j4.15p . we rather hesitate to state this even as a formal asymptotics; further 
mathematical research is necessary and is desirable. 

This Leray's blow-up self-similarity led him to a conjecture on existence of a blow-up similarity 
solution as i — > T^ , and on existence of a self-similar extension for t > T (with (T — t) i— > (t — T) in 
(|5.2p ): see some extra details and a discussion in ^8j § 2]. 
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v^ = Av — \y ■ Vv — 1 V — 




P(v-V)v in M^ X M+, 



where T > is the assumed finite blow-up time of the vector field u(x, t) (used to be 
always T = beforehand). This leads to a different autonomous rescaled equation: 

(5.3) 

where, unlike (15. ip . the box is empty: no exponentially decaying factor therein! 

This is the main principal difference between both key regularity problems: in the 
rescaled equation (15. 3p . the convection term is not accompanied by an exponentially 
small time- factor as in (15. ip . Therefore, it can play a vital role in creating a possible 
L°°-blow-up singularity. It turned out that a self-similar blow-up of Leray's type (15. 2p . 
with a nontrivial profile v = v(y) solving the stationary equation (15. Sp . does not exist: 

(5.4) Av-iy- Vv- iv-P(v- V)v = in R^, v e L'^{M.^) =^ v = 0; 

see [HHl EH EHl [30]. Therefore, it seems that a leading idea how to create a blow-up 
singularity for the NSEs is to deal with the so-called Type II blow-up solutions, i.e., those 
which violate a uniform similarity estimate for Type I solutions: 

(5.5) Type I: for some constant C > 0, |u(x, t)\ < ^^ in R^ x (0, T), 

(5.6) i.e.. Type II blow-up: lim sup \/T — t sup |u(x,t)| = -|-cxd. 

In the rescaled variables (15. 2p . Type II blow-up means looking for a global solution of 
(15. 3p that "blows up" as r — )■ +oo. This problem is open, though some formal scenarios of 
such a Type II blow-up in the NSEs have been discussed for a while; see, e.g., references 
and discussions in [18] . 

Overall, we expect that the present study of simpler boundary point regularity governed 
by an exponentially perturbed rescaled equation (15. ip is an inevitable and important 
step towards solving the main open interior point regularity blow-up problem for (15. 3p . 
Indeed, for (15. ip . it turned out that the asymptotic behaviour near the vertex was spatially 
governed by the simplest first solenoidal Hermite polynomial (actually, a constant). While, 
in order to understand blow-up in the "MPPCI equation" (15. 3p . it seems that an essential 
involvement of all ot/ier solenoidal Hermite polynomials {vo(?/)} as eigenf unctions of the 
linear rescaled operator B* (see Appendix A) should be detected first, before attacking a 
possible essentially nonlinear structure of a Type-II blow-up singularity for the NSEs. 

Acknowledgements. The authors would like to thank I.V. Kamotski for useful dis- 
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Appendix A: Hermitian spectral theory of the linear rescaled operator B*: 

POINT SPECTRUM AND SOLENOIDAL HeRMITE POLYNOMIALS 

Thus, approaching the characteristic vertex (0,0) in the blow-up manner (j2.ip . one observes 
Hermite's operator B* as the principal hnear part of the rescaled equation (J2.2I) . Writing it the 
corresponding divergent form, 

(A.l) B*v = ^V-(p*Vv), 

where the weight is p*{y) = e 4 > 0, we observe that the actual rescaled evolution is now 
restricted to the weighted L^-space L^, (M^), with the exponentially decaying weight p* (y). Here, 
B* is the ("adjoint") Hermite operator with the point spectrum ^ p. 48] 

(A.2) o-(B*) = {Afc = -|, A; = |/3| =0,1,2, ...} (/? is a multiindex in M^), 

where each A^ has the multiplicity '-^ for A^ = 3, or the binomial number C^,^_-^. 

The corresponding complete and closed set of eigenfunctions $* = {ip%{y)} is composed from 
separable Hermite polynomials. Note another important property of Hermite polynomials: 

(A. 3) VV'fl, any derivative D^'ipl is also an eigenfunction with fe = |/3| — I7I > 0. 

Recall that [4] 

(A. 4) polynomial set $* is complete and closed in L^, (M^). 

Further spectral properties are convenient to demonstrate using the linear operator B, 

(A.5) B = A + iyV + |l in Ll(M.^), where p = j^, 

which is adjoint to B* in the dual L^-metric. It has the same point spectrum and the correspond- 
ing eigenfunctions are multiple of the same Hermite polynomials according to the well-known 
generating formula: 

(A.6) ^l;^(y) = t^D^F{y)^rfi{y)Fiy), where F(y) = ^ e'l^'l'/'^ 
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is the rescaled kernel of the fundamental solutions of Dj — A in W^ x ]R_|_. Then, the bi- 
orthonormality holds: 

(A.7) {'il;p,tl;j) =5/3^ for any /3, 7. 

Indeed, this dual metric can be also treated as that in L^*(]R^) for the self-adjoint case, but we 
prefer to keep " L^-dual" notations (for using also in non-symmetric Burnett cases; cf. ()1.7p ). 

Obviously, one needs to consider eigenfunction expansions in the solenoidal restriction 
(A.8) Lj, (M^) = Ll* {R^f n {div v = 0}. 

Indeed, among the polynomials $* = {tp%} there are many that well-suit the solenoidal fields. 
Namely, introducing the eigenspaces 

<5l = Span{i^l, |/3| = k}, k> 1, 
in view of ()A.3p . div plays a role of a "shift operator" in the sense that 
(A.9) div : cDf ^ ^l_^_ 

We next define the corresponding solenoidal eigenspaces as follows: 
(A.IO) Sl = {^* = [vl,v^,v*^f : divv* = 0, < G $*}, where dim 5* = A;(A; + 2); 

see [241 l25l [26] and further references therein. 

Actually, the paper [24] deals with global asymptotics of NSEs solutions as t — >• -|-oo, where 
the adjoint operator B in (jA.SP occurs. Since B is self-adjoint in Lp(M^), many results from 
[25| Append. A] can be applied to B*. For a full collection, see [5l [6] for further asymptotic 
expansions and self-similar solutions. In particular, this made it possible to construct therein 
fast decaying solutions of the NSEs on each ID stable manifolds with the asymptotic behaviouqj 

(A.ll) u^(x,t) ~ t^'="2 v^(^) + ... as t -> 00, where v/3 = v^F G 5fc 

are solenoidal eigenfunctions of B. Namely, taking 

(A.12) ^=[vi,V2,V3f eSk, Vie^k = Span{i;^=^-^D^F{y),\^\ = k}, 

where F stands for the rescaled Gaussian in (|A.6p . we have that 

(A.13) divv = {vi)y, + {V2)y, + {vs)y, = div (v*F) = (divv*)F -^y■^r*F. 

This establishes a one-to-one correspondence between solenoidal eigenfunction classes S^ in 
(jXTOl) for B* and Sk in (jXTT]) for B; see (|XT4l) - (jXT6|) below for the first eigenfunctions 
v/3 = VoF. Therefore, dim Sk = k{k + 2), etc.; see details and rather involved proofs of the 
asymptotics (]A.lip for /c = 1 and 2 in [23]. 

In particular, those solenoidal Hermite polynomial eigenfunctions of B* can be chosen as 
follows [25, p. 2166-69] (the choice is obviously not unique, normalization constants are omitted): 

(A. 14) Ao = : Vg = [1, 1, 1] = e (the first solenoidal Hermite polynomial), 



We present here only the first term of expansion; as usual in dynamical system theory, other terms in 
the case of "resonance" can contain In i- factors {q.v. [1; for a typical PDE application); this phenomenon 
was shown to exist for the NSEs in M^ ^, p_ 236]. 
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(A.15) 



Ai 



'11 



(A.16) 



'21 



'24 





y2 . 

-vl- 

yiy2 

-yivz 



-yiys 
yiy2 . 
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2/3 



-yi. 



'13 



2/1' 



'22 



'25 



2/12/2 

-2/f 
-2/22/3 



2/22/3 



-2/22/1. 



2/i 



-2/2 
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'23 



(dirndl* = 3); 
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-2/22/3 

- 2/? - 2/i. 
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" -2/22/3 " 




2/12/2 




"2/i - 2/i" 


2/22/3 


, V^7 = 


2/1-2/? 


. V^8 = 


-2/12/2 


L2/? - 2/iJ 




. -2/22/3 . 




. 2/12/3 . 



(dimcS2 



etc. 



We need the following final conclusion. By (|A.4p . the set of vectors <^*^ is complete and closed 
iil3ip.(M3)3,sothat 



Ef 



(/3)^/3V^' V 



V* E $f , fe = |/3| > 0, 



=/3V^ = [4(^)^^1'4(^)^^2>4(^)^^3]^- 



(A.17) VvGL2. 

where we use the vector notation 

(A.18) c^(r) = [4(r),4(r),c3(r)f 

It then follows from (jXTjl - lfXH]) that 

(A. 19) polynomial set $* = <&* n {divv = 0} is complete and closed in L^, 

In what follows, we always assume that we deal with "solenoidal" asymptotics involving eigen- 
functions as in (jA.lOp . 

Appendix B: Vertex regularity for Burnett equations 

B.l. Burnett equations in a hierarchy of hydrodynamic models. The Burnett equa- 
tions ()1.7p appear as the second approximation (the NSEs (jl.ip being the first one) of the 
corresponding kinetic equations on the basis of Grad's method in Chapman-Enskog expansions 
for hydrodynamics. Namely, Grad's method applied to kinetic equations, by expanding the 
kernel of the integral operators involved in terms of those with pointwise supports, yields, in 
addition to the classic operators of the Euler equations, other viscosity parts as follows: 

00 
Au = ut + {u-V)u=J2 e2"+^A"(;UnAu) + ... = e(/uoAu + eViA^u + ...) + ... , 

n=0 

where e > is essentially the Knudsen number Kn; see details in Rosenau's regularization 
approach, [62j. In a full model, truncating such series at n = leads to the Navier-Stokes 
equations (jl.ip (with //q > 0), while n = 1 is associated with the Burnett equations (jl.7p . 
Note that Burnett-type equations, with a small parameter appeared as higher-order viscosity 



Note a standard result of functional analysis: polynomials are complete in any weighted L^-space 
with an exponentially decaying weight; see the analyticity argument in Kolmogorov-Fomin j32] p. 431]. 
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approximations of the Navier-Stokes equations, is an effective tool for proving existence of their 
weak ("turbulent" in Leray's sense) solutions; see Lions' monograph [441 § 6, Ch. 1]. Note that 
the "Problem on blow-up/non-blow-up for Burnett equations (|1.7p at an interior point" starts 
from dimensions N = 7: for A'^ < 6, there exists a unique global smooth L^-solution, |2Tl § 6]. 

The final finite-dimensional dynamical system for first Fourier coefficients of solutions to 
(jl.7p is derived similarly, but the regularity conclusions are shown to be more difficult and even 
rather obscure. The necessary spectral properties and vector solenoidal generalized Hermite 
polynomials as eigenfunctions of the rescaled non-self adjoint operator B* are introduced below 
in Appendix C. 

More briefly, we now list below main steps of the regularity analysis for (|1.7p . 

B.2. First blow-up scaling: an exponentially perturbed parabolic equation. The first 
blow-up scaling in (11. 7p is now 



(B.l) u(x,t)=v(y,r), y = jz^, r = - In(-t) : (-1,0) ^ M+, 

which yields the following rescaled equation: 

(B.2) v^ = -e"xVp + B*v-e~x (v V)v in Qo, where B* = -A^ - i y • V 

is the adjoint operator (jC.ip for m = 2 with good spectral properties given in Appendix C. 

B.3. Backward paraboloid. Here, dQo is defined as follows: 

(B.3) S{t) = ^Q,r^{t}■. qo{x)^{Y.t,a,xf)-- = {-t)\ip{T) {Y.< = 1\ 

with the same slow growing functions '^{t) as in ()2.3p . The rescaled equation ()B.2p is then set 
in an expanding domain, 

(B.4) 5(t) = 5Qo n M : q^{y) = [Ylti^ivfY = ^i-^) ^ +'^ as r ^ +oo. 

B.4. The Cauchy problem setting. Extending to the Cauchy problem by using the variables 
(I2.12p . we use Green's second formula: for any x ^ C'o"; 

(B-5) 4o(.)<.M} (-^'X - X^'-)dy = hit) ( Axf^ - v^ - Avg + x^) ds. 

In view of the Dirichlet boundary conditions in (jl.7p . in the sense of distributions, ()B.5P reads: 

(B.6) A^v = A^v// + (|,Av) <5^(^) + £(Av5^(^)), 

where densities of a single- and a double-layer potential now depend on the Laplacian Av instead 
of V in ()2.12p . Using the same harmonic pressure extension, we obtain 

(B.7) V. = B*v - e-^ P (v • V)v - P (|^Av) 5^^^^ - P |i(Av <5^(,)). 

As for the NSEs with m = 1, this problem is always locally well-posed, and is guaranteed to be 
globally well posed either for A^ < 6 or for any sufficiently smooth initial data for N >7. 
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B.5. Boundary layer and a perturbed rescaled equation. Close to the lateral boundary 
of Qo) the next rescaled variables are 

(B.8) ^^W) ^"^ v(y,r) = w(z,t). 

This makes the corresponding rescaled paraboloid ()B.4p fixed: 
(B.9) 5: {Y^^,a,zf)-^=l. 

The rescaled vector field w satisfies a perturbed equation: 

w, = -^A2w-iz.V,w+^z.V,w-ie-^P(wV,)w 

The BL-variables for a fixed point ^o ^ 5* on the boundary ()3.2p are 

4 4 

(B.ll) i = ip^i{T){zQ- z), 993(r)dr = ds, and w(2;,t) = p(s)g(C, s). 

In this boundary layer, we are looking for a generic pattern of the behaviour described by (|B.10p 
on compact subsets near the lateral boundary, satisfying 

(B.12) 1^1 =o((/?"i(r)) ^0 =^ |z-zo| =o((/5-i(r)) ^0 as r ^ +00. 

Substituting (IBTT]) into the PDE (iBlOl) yields 

gs = Ag - i ^ ^ V^g - f (zo - ^) • V^g 
(B.13) -I J^ ^ V,g - f g + I e-^ F(g • V,)g - i ^^-^-t P (|,Av) 5^ 

_i^iV-4-|p|.(Av5^), where Ag = -A^g + i zq • V^g. 

Again, in ()B.13p . we observe a perturbed linear uniformly parabolic equation. As usual, here 
one needs to check that all the perturbation terms are asymptotically small as s (or r) — )• +00, 
relative to the stationary autonomous operator A. This is done similarly to m = 1 above. 
Overall, the BL representation (jB.lip and (j3.5p imply that ()3.8p holds. Moreover, asymptotically, 
the limit problem becomes one-dimensional, depending on the space variable (j3.9p . 

We again pose the same asymptotic behaviour (13. 8p at infinity. According to the scaling 
(jB.lip . let us fix a uniformly bounded rescaled orbit {g(s), s > 0]lj. Then, by parabolic theory 
[U [To], we can again pass to the limit in (IB.13P along a subsequence {sfc} — )■ +00, removing 
small perturbations. Therefore, uniformly on compact subsets defined in (IB.12p . as /c — )• 00, 

(B.14) g(sfc + s)^h(s), where h, = Ah, h = f^ = at ^ = 0, hi\^=+^ = l. 

The liinit equation obtained from ()B.13p . 

(B.15) hs = Ah = -\rir,'n + 3 ^»? 

is again a standard linear parabolic PDE in M^ x M^., with a non self-adjoint operator A, so 
()B.15P is not a gradient system in L?. We then need to show that, in an appropriate weighted 



As usual, the scaled function p{s) remains unknown and to be determined by matching with the 
inner region behaviour. 
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-L^-space if necessary and under the hypothesis (j3.8p . the stabihzation holds, i.e., the w-hmit set 
of the orbit {h(s)}s>o consists of a single equilibrium: as s — )■ +00, 

{h(^,s) —?- goiO, where Ago = for r/ G M, 
go = g'o = for r/ = 0, gii+oo) = 1. 
This gives the unique solution of (|B.16p (see [T71 § 7] and [201 § 5]): for, e.g., zq = [1, 1, 1]-^, 
(B.17) 5^(^) = i_e-W^[cos(^) + ^sin(^)], j = 1, 2, ... , TV. 

It turns out that the limit problem (JB.ISP possesses a number of strong gradient and contrac- 
tivity properties. Namely setting by linearization 

(B.18) h(s) = go + w(s) =^ Ws = Aw = -w^^^^ + j w^, w = w^ = at r/ = 0, 

we arrive at the following (cf. Proposition 13. II for m = 1): 

Proposition B.l. (i) (jB.lSP is a gradient system in L"^ , and 

(ii) In the given class of solutions, the oj-limit set J^o of (]B.18P consists of the origin only and 
is uniformly stable. 

Proof, (i) One can see that (jB.lSp admits a monotone Lyapunov function obtained by multiply- 
ing by w^^ in L^: 

(B.19) \^J\{w%Y = -!{{w%,,Y<Q. 

Hence, (ii) also follows. D 

Thus, quite similar to the second-order case, under given assumptions, we can pass to the limit 
s — )• -|-oo along any sequence in the perturbed gradient system (jB.lSp . Then, again similarly to 
m = 1, the uniform stability of the stationary point go in the limit autonomous system (jB.lSP in 
a suitable metric guarantees that the asymptotically small perturbations do not affect the omega- 
limit set; see |23t Ch. 1]. However, at this moment, we cannot avoid the following convention, 
which for m = 2 is much more key than for m = 1. Actually, the convergence ()B.14p and ()B.16P 
for the perturbed dynamical system (IB.ISP should be considered as the main Hypothesis, which 
characterizes the class of generic patterns under consideration, and then the normalization (|3.8p 
is its partial consequence. For bi-harmonic flows, a more clear characterization of this class of 
generic patterns is very difficult. It seems that a correct language of doing this (in fact, for both 
cases m = 1 and tti > 2) is to reinforce the corresponding centre subspace behaviour as in (j4.7p . 

Finally, we summarize these conclusions as follows: 

Proposition B.2. Under the given hypothesis and conditions, the problem (jB.lSp admits a 
family of solutions (called generic) satisfying (|B.16p . 

Such a definition of generic patterns looks rather non-constructive, which is unavoidable for 
such higher-order nonlocal PDEs. However, ()B.16P is expected to occur for "almost all" solutions. 

Thus, we stop further discussions concerning the passage to the limit s — )• -|-cxd in (]B.13p . 
which, as we have shown, under the given hypotheses on the asymptotic smallness of pertur- 
bations available, reduces to a linear stability analysis of the nontrivial equilibrium go of the 

27 



linear rescaled operator A in (|B.13p . This has been resolved for a class of generic solutions. 
More generally, we have to deal with solutions of ()B.13p from the stable subset Wo of go within 
the prescribed perturbed equations (jB.lSp . A clear, constructive, and full identification of Wo 
is not possible for such higher-order nonlocal perturbed parabolic equations. 

We summarize the conclusions as follows: in what follows, under the given hypothesis and 
conditions, we will deal with a family of solutions Wq of (JB.ISP {called generic), for which 
(lBl6]) holds. 

B.6. Inner region analysis. In the Inner Region, the original rescaled problem ()B.7p occurs. 
For any extended solution orbit (j2.1ip uniformly bounded in L'i* (M) , we use the eigenfunction 
expansion ()4.3p via the generalized solenoidal Hermite polynomials (IC.2ip . Substituting (14. Sp 
into ()B.7p and using the orthonormality property ()C.20p yield a dynamical system: for any 
multiindex /3, with \/3\ > 0, 

(B.20) {c^ = X^cp - (P (I^Av) 5^(,), v;,) - (P|^(Av5^(^)), v^) - e'X (P(v • V)v, v^), 

where A^ = —-^ by ()C.14p . so that A/3 < for any |/3| > 1. As for m = 1, one then needs to 
concentrate on the first Fourier generic patterns associated with the centre subspace for B* (cf. 
(fOMP for iV = 3) 

(B.21) k = 0: Ao = and v*o{y) = e = [1,1, ...,lf , Mv) = F{y)e. 

This reflects another characterization of our class of generic patterns. The equation for co(t) 
then takes the form: 

(B.22) co = - / (I^Av)vods- / Av|^vods-e-^ / P(v • V)v vq dy. 

Sir) Sir) M.N 

Note that first two terms on the right-hand side in ()B.22p have a pure solenoidal parabolic 
(bi-harmonic) nature, while the only Navier-Stokes influence is presented by the last nonlinear 
term with an exponentially decaying factor. 

Using next the boundary behaviour ()B.16P with the ID profile (JB.ITP for r >> 1: in the 
rescaled sense, on the given compact subsets, (|3.17p holds, with go given by (jB.lTp . where rj 
stands for the rescaled distance: 

(B.23) r? = v9^(r)dist{y,5(r)}. 

Similar to ()3.18p . such a BL-asymptotics is "almost" solenoidal for r ^ 1, i.e., perturbed by an 
exponentially small factor at any distance (^0 > from the boundary. By the matching of both 
Regions for such generic patterns, (j4.8p has to remain valid. 

Performing, similar to (j4.10p - (|4.19p . proper estimating of all the three terms on the right- 
hand side of ()B.22p (the last, the Navier-Stokes one can be again estimated rather roughly) 
yields the following dynamical system for the first Fourier coefficients: 

(B.24) |co = -7iiCo(/7^vo((/9) - 7i2CoV'^" ^ Vo((/7) - 72 e"^ (cq • e)co 93^+^ vo{ip) + ... , 

where, again, the first two terms are purely "parabolic". Similar to (j4.20p . we replace surface 
integrals by some "average" values, actually assuming the radial dependence on \y\ with the 
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rescaled surface 5(r) : {\y\ = (/?(r)}, as in (j4.13p . As above, we do not guarantee that the 
niultipher 99^+^ in the last term in (]B.24p is any optimal one (since the whole term will be 
shown to be negligible anyway). 

Finally, using the expansion of the rescaled kernel given in (jC.lip and keeping the leading 
term only, the asymptotic dynamical system reads 
/■ 



(B.25) 

-72 e 



^ (Co •e)co(^^+-^-^»(r)e-'^«^*'"W[C7isin(6o(/'^(r)) + C2cos(6ovpt(r))], 



where, as usual, 71,2 £ 1^^ are some constant vectors. 

We then arrive at a typical and simple ODE regularity criterion: under the given hy- 
potheses and conventions of our asymptotic analysis, the vertex (0, 0) is regular in the class of 
generic solutions iff any solution of the non- autonomous 3D dynamical system ()B.25P vanishes 
as T ^ +00, i.e., is globally asymptotically stable for ()B.25p . 

B.7. Two regularity conclusions. 1. The first regularity conclusion is straightforward: in 
the absence of the convection term, i.e., for the linear fourth- order Stokes problem, 

(B.26) ut = -Vp - A^u, divu = in Qq, 

with the Dirichlet boundary conditions as in (11. 7p . the vertex (0,0) is regular provided that the 
following integral diverges to —00 (cf. (j2.9p ): 

(B.27) 



-t-00 
^71 / 99^+ 3 -^« (s) e-'^«'^'^'(^) [Ci sin(6oV?3 (s)) + C2 cos(6oV'^ («))] ds 



Since the function under the integrals is strongly oscillatory and, in general, is of changing sign 
for s 3> 1, ()B.27p may require a special procedure of "oscillatory cut-ofF' of the given (^(r) 
to delete a possible positive part of the divergent integral; see [iTl § 7]. In other words, to 
get the regularity conclusion ()B.27p . the behaviour of v?(t) as r — ?• +c« must be very carefully 
adjusted with the nonmonotone and oscillatory behaviour of the rescaled kernel F((^(t)) of the 
fundamental solution ()C.4p of the bi-harmonic operator. 

If ()B.27p fails, then the vertex is not regular. One can see that the transition from a possible 
regularity (after an oscillatory cut-off) to the guaranteed irregularity occurs at the following 
"critical" paraboloid with 
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3 11 3 „ ,3 



(B.28) doipS{T) = lnT =^ lp^{t) = d^^ (in t)^ = 2- 3~^ (In t)* as r ^ +00. 

This is the fourth-order analogy of Petrovskii's function ()2.8p . so that the constant therein, 
(B.29) a = 2X 3" ' 
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is optimal (similar to the "2" in (12. Sp ): replacing it by any larger one C* -|- e, with an e > 0, 
guarantees convergence in (|B.27p and hence the irregularity of the vertex (0, 0). Indeed, for such 
99(r), the integral in ()B.27p simply converges, i.e., on almost all such centre subspace orbits, 

(B.30) co(r) 7^ as r ^ -Foo. 
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2. Concerning the full nonlinear Burnett problem (jl.7p . we arrive at the same conclusion as 
for the NSEs. Namely, it follows by balancing two terms on the right-hand side of ()B.25p . that 
the last nonlinear one may be leading provided that 

(B.31) |co(t)| >e^ ^3^:;^ ^+00 as r ^ +00, 

which never happens in the case of the "linear" regularity described by ()B.27p . Hence, the 
convective term cannot change the vertex regularity, thus leading to the same regularity criterion. 

B.8. Nonexistence of similarity blo^v-up for Burnett equations: Type II singularities 
are also needed. For ()1.7p . Leray's-type blow-up scaling ()5.2p takes the form 

In(T-t), 



(B 


.32) 




u(x,t) = (T-t)-4v(y,r 
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v^ = B*v - f V - 
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'V)v 


in 


R^ 


X 



Here B* is the linear rescaled operator (jC.lh with the known point spectrum and eigenfunctions 
being generalized Hermite polynomials. 

As for the NSEs (jl.ip (see Section [5]), the first question is whether a nontrivial Type I self- 
similar blow-up exists, i.e., whether a nontrivial stationary solution v = v(y) of (1B.33P exists: 



(B.34) -A^v-iyVv- |v-P(v V)v = in R^, v G L^ 

It is curious that a negative answer (i.e., similar for the NSEs in Section[5]) can be obtained rather 
convincingly just by a local asymptotic analysis of the elliptic equation (|B.34p . As happens in 
practically all blow-up problems for reaction-diffusion and other nonlinear PDEs (see examples 
in, e.g., [iHl [23l [63] ) , a "generic" behaviour of its solutions as z = \y\ — )• +00 is governed by the 
leading lower-order linear terms, i.e., in the radial representation, this ineans that, for z ^ 1, 

(B.35) -izv^ - |v + ... = =^ v(2;) ~ ^ as z ^ +00, 

where C € M^ is a constant vector. Of course, (JB.35P is just a rough radial estimate, so an extra 
"angular separation" is necessary to produce all asymptotics like that at infinity. However, 
()B.35P is sufficient for a key negative conclusion: via the local behaviour (|B.35p . for any C 7^ 0, 

(B.36) ^ e L2({|y| > 1}) iff N<6. 

In other words, in the "blow-up case"L3 N > 7, blow-up cannot be of a self-similar (Type I) 
form (|B.32p with a nontrivial asymptotics ()B.35p . 

Surely, this is not a proof of such a nonexistence, since the special single case C = in (IB.35P 
has not been ruled out. Indeed, formally, it can happen that, for C = 0, the similarity profile 
v(y) solving ()B.34p inay reach an exponential decay at infinity (on derivation, see Appendix C) 

(B.37) v(y) ~Ci|iye-"ol2'l*''', where oq = 3-2-1 and Ci G M^. 



l^For iV < 6, solutions of dUT]) do not blow-up in i°°(]R^), 

30 



Example: a diversion to blo^v-up in a related semilinear bi-harmonic flo^v. As is 

known from similarity blow-up in semilinear bi-harmonic equations such as 

(B.38) ut = -A^u+\u\P''^u in M^ x M+, where p>l, 

such a behaviour is highly unlikely. To explain this, consider its self-similar blowing-up solutions 

(B.39) u{x, t) = {T- t)-^ v{y), y = j^^ =^ -A^v -ly.Vv-^v + \v\p~^v = 0. 

Checking the asymptotic behaviour as y — )• -|-oo, we again obtain the generic algebraic decay 
similar to ()B.35p . which is governed by two linear terms: as z = \y\ — )• -|-oo, 

(B.40) -lzv'-^v + ... = =^ v{y)^C\y\-^K 

Similarly, for C = 0, the behaviour gets exponentially decaying (cf. ()B.37P ): 

(B.41) i;(y) ~ Ci|y|V"ol2^l'''', where ,^ = -|(iV-^), 

with the same oq as in ()B.37p . See [191 § 2.3] and Appendix C for a derivation of such two- 
scale WKBJ-asymptotics ()B.4ip and ()B.37p . Therefore, the asymptotic bundle of exponentially 
decaying solutions (lB.4ip of the ODE in ()B.39p contains a unique parameter Ci E M, i.e., it is 
one-dimensional. Thus, this is not enough to "shoot" two symmetry conditions at the origin: 

(B.42) v'{0) = v"'{0) = 0, 

so an extra parameter should be at hand, and this is p. As shown in [191 § 2] by a careful 
numerical analysis of the ODE in ()B.39P for A^ = 1, there exists a unique value of the exponent 

(B.43) p = ps = lAO..., 

for which (IB.39P admits a solution with the exponential decay ()B.4ip . with some Ci ^ 0. Then, 
in the sense of bounded measures in M, for p = ps, 

(B.44) \u{x,t)\'^ ^ D6{x) as t ^ T^ , with the constant Z:* = / |f(y)|^ dy > 0. 

Back to blow-up in Burnett equations. We expect that a similar phenomenon does not 
exist for similarity blow-up in the Burnett equations, i.e., exponentially decaying similarity 
profiles ()B.37p do not exist. Recall that, unlike ()B.38p . the equations (|1.7p and ()B.34p do not 
contain any free parameter (like p in (IB.38P ) , which could allow to get such a solution at least for 
some its values. Of course, (IB.34P is a system of three solenoidal fourth-order semilinear elliptic 
equations, and a definite negative nonexistence conclusion is very difficult to justify rigorousljo. 

Overall, we arrive at the following plausible situation: similar to the NSEs (jl.ip in dimensions 
N > 3 {see discussion in Section\^, blow-up in the Burnett equations (11. 7p in dimensions N >7 



An extra parameter may be "hidden" in a kind of "symmetry group" in tlie R^-geometry admitted 
by these PDEs (anyway, this looks not that convincing). Overall, existence of a pure self-similar blow- 
up for Burnett equations for iV = 7 is a too simple way to settle this new "fourth-order Millennium 
Problem", and (at least one of the) authors would like to rule out such a trivial solution of it. 

31 



cannot be self-similar and requires constructing (or proving their nonexistence) nan- self- similar 
Type II blow-up singularities i. 



Appendix C: Solenoidal Hermitian spectral theory for operator pair {B, B*} 

We describe the necessary spectral properties of the hnear 27Ti,th-order differential operator 
in M^ {m = 2 for the Burnett equations ()1.7p ) 

(C.l) B* = i-ir+^A^-^y.Vy, 

and of its L^-adjoint B given by 

(C.2) B = (-ir+iA- + ^y.V, + ^/. 

As we have seen, for m = 1, (IC.ip and (IC.2P are classic Hermite self-adjoint operators with 
completely known spectral properties, [H p. 48]. For any m > 2, both operators (jC.ip and 
(IC.2p . though looking very similar to those for m = 1, are not symmetric and do not admit a 
self-adjoint extension, so we follow [8j in presenting spectral theory. 

C.l. Fundamental solution, rescaled kernel, and first estimates. The fundamental so- 
lution b{x,t) of the linear poly-harmonic parabolic equation 

(C.3) ut = -{-A^u in M^ x M+ 

takes the standard similarity form 

(C.4) b{x,t)=t-^F{y), y = ^. 

The rescaled kernel F is the unique radial solution of the elliptic equation 

(C.5) BF = -{-ArF + ^yVF + ^F = mR^, with J F = 1. 

For m>2, the rescaled kernel function i^'dy]) is oscillatory as \y\ — )• oo and satisfies [10^ 113] 

(C.6) \F{y)\ < D e-"'«l?'l° in M^, where a = ^^ G (1,2), 

for some positive constants D and do depending on m and N. 



Here, there occurs a 4th-order Blow-up Problem for the Burnett equations (|1.7p that may be much 
more difficult mathematically than the Millennium Prize Problem for the Navier-Stokes ones (jl.ip . Of 
course, unlike the classic one in the actual M.^, a "non-realistic" dimension N = 7 makes the 4th-order 
Problem less attractive for applications and for a general public, but, mathematically, it can be even more 
fundamental for PDE theory, since represents less understood features and principles of interaction of a 
higher-order viscosity-diffusion operator with a nonlinear convection one gathered in a nonlocal fashion. 
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C.2. Some constants. As we have seen, the rescaled kernel F{y) satisfies (|C.6p . where d^ 
admits an exphcit expression; see below. Such optimal exponential estimates of the funda- 
mental solutions of higher-order parabolic equations are well-known and were first obtained by 
Evgrafov-Postnikov (1970) and Tintarev (1982); see Barbatis 121 Sj for key references. 

As a crucial issue for the boundary point regularity study, we will need a sharper, than given 
by ()C.6p . asymptotic behaviour of the rescaled kernel F{y) as y — )• +oo. To get that, we keep 
four leading terms in (IC.5J1 and obtain, in terms of the radial variable y i— )• \y\ > 0: 

(C.7) (-l)™+i[F(2m)+^A^^(2m-i)^^_^j ^_I.y^/ + ^^ = for y » 1. 

Using standard classic WKBJ asymptotics, we substitute into ()C.7p the function 

(C.8) F{y) = y-^^ e'^y^ + ... as y ^ +00, 

exhibiting two scales. Balancing two leading terms gives the algebraic equation for a and 6q: 

(C.9) (-l)-(aa)2-i = A. and <^o = "^^TJ}-^ > • 

By construction, one needs to get the root a of (|C.9P with the maximal Re a < 0. This yields 
(see e.g., [21 [3] and [22l p. 141]) 

(CIO) a=^[-,m{^^^^)+icos{^^^^)]^-do + ibo (do > 0). 

Finally, this gives the following double-scale asymptotic of the kernel: 

(C.ll) F{y) = y--^" e-'^oy°' ^Ci sin(6oy") + C2 cos(6oy")] + ••• as y = \y\ -^ +00, 

where Ci,2 are real constants, |Ci| + IC2I 7^ 0. In (IC.lip . we present the first two leading terms 
from the TTi-dimensional bundle of exponentially decaying asymptotics. 
In particular, for the Burnett equations (]1.7p in M^, we have 

(C.12) m = 2,N = 3: a = |, do = 3- 2"^, 60 = 35 -2"^, and ^o = |- 

C.3. The discrete real spectrum and eigenfunctions of B. For m > 2, B is considered in 
the weighted space -Lp(M^) with the exponentially growing weight function 

(C.13) p(y) = e-^l^l" > in M^, 

where a G (0, 2do) is a fixed constant. We next introduce a standard Hilbert (a weighted Sobolev) 
space of functions Hp"^(M.'^) with the inner product and the induced norm 

2m 2m 

{v,w)p = J p{y) Yl D^v{y) D^wiy) dy and \\v\\l = J p{y) ^ \D'^v{y)\'^ Ay. 

RJV A:=0 RAT fc=0 

Then Hf^{^^) C ^^(M^) c L2(M^), and B is a bounded linear operator from ij2™(M^) to 



Lp(R^). Key spectral properties of the operator B are as follows [8]: 

Lemma C.l. (i) The spectrum ofB comprises real simple eigenvalues only, 

(C.14) a(B) = {A^ = -^, k = |/3| =0,1,2,...}. 

(ii) The eigenfunctions i^piy) are given by 

(C.15) V/3(y) = ^^ D^Fiy), for any \^\=k> 0. 
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(iii) Eigenfunction subset (|C.14p is complete in L'^{]R.) and in Lil 

(iv) The resolvent (B — XI)^"^ for A o"(B) is a compact integral operator in L?{M. ). 

By Lemma IC-H the centre and stable subspaces of B are given by 

(C.16) E^ = SpanjV'o = F] and E' = Spanj^-/?, |/3| > 0}. 

C.4. Polynomial eigenfunctions of the operator B*. Consider the operator (jC.ip in the 
weighted space L^, (M ), where {■,-)p* and || • \\p* being the inner product and the norm, with 
the "adjoint" exponentially decaying weight function 

(C.17) p*(y)^-^=e-«H">0. 

We ascribe to B* the domain ff^™(]R^), which is dense in L^.(M-'^), and then 

B* : H^p^{R^) ^ Ll,{R^) 

is a bounded linear operator. B is adjoint to B* in the usual sense: denoting by (•, •) the inner 
product in the dual space L^(]R'^), we have 

(C.18) {Bv,w) = {v,B*w) for any t; G F2"^(M^) and u; G ^^^(R^). 

The eigenfunctions of B* take a particularly simple finite polynomial form and are as follows: 

Lemma C.2. (i) cr(B*) = a{B). 

(ii) The eigenfunctions ^'siv) o/B* are generalized Hermite polynomials of degree |/3| given by 



(C.19) reiy) = :4 k + EKf "' U-^y^'y 



VW- 



-i=i j\ ' 



for any (3. 



(iii) Eigenfunction subset (JC.IOP is complete in L^, (M ). 

(iv) B* has a compact resolvent (B* - A/)"^ in L^.(R^) for A cr(B*). 

(v) The bi-orthonormality of the bases {ipii\ and {V'X} holds in the dual L^-metric: 

(C.20) (V'/3,V'7) = <J/37 for any f3, 7. 

Remark on closure. This is an important issue for using eigenfunction expansions of solutions. 
Firstly, in the self-adjoint case m = 1, the sets of eigenfunctions are closed in the corresponding 
spaces, [4] (and we have used this in our previous NSEs study). 

Secondly, for m > 2, one needs some extra details. Namely, using ()C.20p . we can introduce 
the subspaces of eigenfunction expansions and begin with the operator B. We denote by L^ the 
subspace of eigenfunction expansions v = Yl^^pi^P with coefficients c^ = {v,tp*) defined as the 
closure of the finite sums {X]|«|<^f c^V'^} in the norm of L^. Similarly, for the adjoint operator 
B*, we define the subspace L^* C L^*. Note that since the operators are not self-adjoint and 
the eigenfunction subsets are not orthonormal, in general, these subspaces can be different from 
Lp and L^, , and particularly the equality is guaranteed in the self-adjoint case m = 1, a = j. 

Thus, for m > 2, in the above subspaces obtained via a suitable closure, we can apply standard 
eigenfunction expansion techniques as in the classic self-adjoint case m = 1. 
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C.5. Solenoidal Hermite polynomials. The vector solenoidal Hermite polynomials are con- 
structed from (IC.19P in a manner similar to that for m = 1; ci (lA.14p -( lA.16p . Namely, given a 
vector polynomial 

(C.21) v; = [V;,,V;„...,V'^^f, where |/3i| = |/32| = ... = |/3^| = |/3|, 

it gets solenoidal provided that 

N 

(C.22) divv;=E(#Jj/. =0. 

For instance, for the Burnett case m = 2 and A^ = 3, some pairs are (not all linearly indepen- 
dent eigenfunctions are presented, normalization constants are omitted): 

(C.23) Ao = 0: vS = [l,l,lf, 

(C.24) Xi = -l-^u = [y2,-y3,y2f, v^2 = [2/3,2/3, -yi]'^, '^h = [-y2,yi,yiV, 

(C.25) A2 = -^: ^2i = [-yi-y3^yw2,yiy3f, ^22 = [yiy2,-y2-y3^ y^y-if, etc. 

(C.26) A3 = -|: v*i = [y|,y|,y3], v*2 = [yiyi,y2yi, -2/3(2/? + y|)], etc. 

(C.27) A4 = -i: ^^l, = [yt + A\,yl + 4l,yf + Ay]^, ^l2 = [yiyly2yf,-y3{yf + yl)f , etc. 

As in the self-adjoint case m = 1, some technical efforts are necessary toward complete- 
ness/closure of generalized solenoidal Hermite polynomials in suitable spaces. We omit details. 
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